Density of critical clusters in strips of strongly disordered systems 
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We consider two models with disorder dominated critical points and study the distribution of 
clusters which are confined in strips and touch one or both boundaries. For the classical random 
bond Potts model in the large-g limit we study optimal Fortuin-Kasteleyn clusters by combinato- 
rial optimization algorithm. For the random transverse-field Ising chain clusters are defined and 
calculated through the strong disorder renormalization group method. The numerically calculated 
density profiles close to the boundaries are shown to follow scaling predictions. For the random 
bond Potts model we have obtained accurate numerical estimates for the critical exponents and 
demonstrated that the density profiles are well described by conformal formulae. 
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I. INTRODUCTION 

In a critical system the correlation length is divergent 
and correlated domains appear in all length scales. This 
phenomenon is seen for percolation 1 - where the domains 
are connected clusters. In discrete spin models, such as 
in the Ising and the Potts models, domains of correlated 
spins can be identified in different ways. One possibility 
is to use geometrical clusters^ (also called Ising or Potts 
clusters) which are domains of parallel spins. In two di- 
mensions (2d) geometrical clusters percolate the sample 
at the critical temperature and their fractal dimension 
can be obtained through conformal invarianco 3 -. This 
value is generally different from the fractal dimension of 
Fortuin-Kasteleyn (FK) clusters^ which are represented 
by graphs of the high-temperature expansion. From a 
geometrical cluster the FK cluster is obtained by remov- 
ing bonds by a probability, 1 — p = e~ K ", K c being the 
critical value of the coupling. The fractal dimension of a 
FK cluster is directly related to the scaling dimension of 
the magnetization. 

In a finite geometry, such as inside strips or squares, 
one is interested in the spanning probability and dif- 
ferent crossing problems of the critical clusters. For 
2d percolation many exact and numerical results have 
been obtained in this fi e ldii£iL&MJiIiii2iIii. Another in- 
teresting problem is the density of clusters in restricted 
geometries 1 ^, which is defined by the fraction of sam- 
ples for which a given point belongs to a cluster with 
some prescribed property, such as touching the edges of 
infinite and half-infinite strips, squares, etc. This lat- 
ter problem is analogous to the calculation of order pa- 
rameter profiles in restricted geometries, which has been 
intensively studied through conformal invariance and nu- 
merical methods 1 ^^^^^^^^^ 3 .^^^^^ 2 ^. 

Correlated clusters can be defined also in models in 
the presence of quenched disorder. In a random fixed 
point one generally considers such quantities, which are 
averaged first over thermal fluctuations and afterwards 
over quenched disorder. In isotropic random systems 



conformal symmetry is expected to hold at the critical 
point so that average operator profiles and average clus- 
ter densities are expected to be invariant under conformal 
transformations. Among disordered systems an interest- 
ing class is represented by such models in which the tran- 
sition in the pure version is of first-order, but in the disor- 
dered version the transition softens into second order—. 
This type of random fixed point can be found, among 
others, in the two-dimensional random bond Potts model 
(RBPM) for q > 4, q being the number of states^L^ 2 .. 

If the distribution of the disorder is not isotropic, e.g. 
it has a layered structure, then the scaling behavior of the 
disordered system can be anisotropic, which is manifested 
in the fact that the critical clusters have an elongated 
shape. This means that the characteristic sizes of the 
clusters parallel, and perpendicular, to the lay- 
ers are generally related as £ii ~ with an anisotropy 
exponent, z / 1. These essentially anisotropic models 
arc not conformally invariant. A well known example 
in this class is the McCoy- Wu model^ 3 -, which is a two- 
dimensional Ising model with layered disorder. Study of 
this system, as well as its one-dimensional quantum ver- 
sion the random transverse-field Ising chain (RTFIC) has 
shown^ 4 - that the critical behavior is controlled by a so 
called infinite disorder fixed point (IDFP), in which scal- 
ing is strongly anisotropic^. The characteristic lengths 
are related as ln£|| ~ so that the anisotropy exponent 
is formally infinite. The same IDFP is found to con- 
trol the critical behavior of the randomly layered g-state 
Potts model 3 ^, as well as for strong enough layered dis- 
order the critical behavior of percolation 3 ^ and directed 
percolation^. Operator profiles in the RTFIC have been 
studied numerically 2 ^ and the obtained data could be 
well fitted by curves which are obtained by analogy of 
the conformal results. 

In this paper we study the density of critical clus- 
ters in two problems in which the critical properties are 
dominated by strong disorder effects. The first model 
is the two-dimensional RBPM in the large-g limit. In 
this model for a given realization of disorder the high- 
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temperature series expansion is dominated by a single 
graph^ the so called optimal diagram, thus thermal 
fluctuations are indeed negligible. This optimal dia- 
gram is calculated for each finite sample by a combi- 
natorial optimization algorithm 4 ^. Clusters in the opti- 
mal diagram are isotropic and density of clusters is ob- 
tained through averaging over disorder realizations. The 
second model we consider is the RTFIC, i.e. a ran- 
dom quantum model which is related to the classical 
McCoy- Wu model, in which the Fortuin-Kasteleyn clus- 
ters are strongly anisotropic. In the RTFIC clusters of 
correlated spins can be defined and calculated by the 
so called strong disorder renormalization group (SDRG) 
method^. During renormalization the system is trans- 
formed into a set of effective spin clusters and for a finite 
system with a given realization of the disorder one ob- 
tains the final cluster, which contains the mostly corre- 
lated sites. The fractal dimension of the final cluster at 
the critical point is directly related to the scaling dimen- 
sion of the magnetization of the RTFIC. Here we calcu- 
late the density of these final clusters which are confined 
in a (one-dimensional) strip. 

The two models we study in this paper are expected 
to be closely related, as far as their critical proper- 
ties are concerned. Based on numerical and analytical 
studie s 41 ' 42 the scaling dimension of the magnetization, 
Xb, and that of the surface magnetization, x s , are conjec- 
tured to be the same for both systems and given by^: 



3-^5 _ 1 

Xb — ^ > x s — — 



(1) 



On the other hand the correlation length exponents are 
related by a factor of two: it is v = 2 for the RTFIC and 
v = 1 for the RBPM. Here we are interested in a possible 
analogy in terms of the densities of the critical cluster. 

The structure of the paper is the following. SecHIl is 
devoted to the RBPM. Here we define the model, outline 
the calculation of the optimal diagram and then ana- 
lyze the statistics of the distribution of the clusters. The 
numerically calculated densities are then compared with 
formulae which are obtained by modifying conformal re- 
sults for percolation. In Scc lIIII we define the RTFIC, 
recapitulate the essence of the SDRG method and then 
numerically calculate final clusters at the critical point. 
The numerically calculated densities are compared with 
analytical formulae in this case, too. The paper is closed 
with a discussion. 



II. OPTIMAL CLUSTERS IN THE RBPM 
A. Model 

The g-state Potts model 4 ^ is defined by the Hamilto- 



h = ~y] Jij5(<Ji,<jj) 



(2) 



in terms of the Potts-spin variables, = 0, 1, ■ • ■ , q — 
1, at site i. The summation runs over all edges of a 
lattice € E, and in our study the couplings, > 
0, are independent and identically distributed random 
numbers. To write the partition sum of the system it is 
convenient to use the random cluster representation 4 -: 



z = E « C{G) II b" Jii i] 



(3) 
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where (3 = 1/ (ksT In q) , the sum runs over all subset 
of bonds, G C E and c(G) stands for the number of 
connected components of G. In the following we restrict 
ourselves to the square lattice, in which case the phase 
transition in the non-random model is of second order 
(first order) for q < 4 (q > 4)^, but for random couplings 
the phase transition softens to second order for any value 
of g 45 i 46 . For conceptional simplicity we consider the 
large-g limit, where q@ Ji i 3> 1, and the partition function 
can be written as 



GCE 



<t>(G)=c(G)+pJ2 J H ( 4 ) 
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which is dominated by the largest term, <j>* = maxg 4>{G). 
Consequently for a given realization of disorder the ther- 
mal fluctuations play a completely negligible role and the 
critical properties of the system arc dominated by disor- 
der effects. The optimal diagram of the RBPM plays 
a completely analogous role as the geometrical clusters 
in percolation theory. For example at the critical point 
there is a giant cluster in the optimal diagram the fractal 
dimension of which, df, is related to the scaling dimen- 
sion of the (average) magnetization as d = dt+Xb, where 
d = 2 is the dimension of the system. One can also study 
other questions, such as distribution of the mass of the 
connected clusters, spanning probability, surface scaling 
exponent, etc. Here we are going to investigate the den- 
sity of clusters in strip geometry. 

During our study we use a bimodal form of the dis- 
order, when the reduced couplings, Kij = /3Jij take two 
values: K\ = K — A and K2 = K+A with equal probabil- 
ity. Generally we study the critical point of the system 
which is located at K = K c = 1/2^ independently of 
the value ofO < A < 1/2. Note that the pure system 
is obtained for A = 0, whereas for A = 1/2, when just 
the strong bonds are present in the system we have the 
traditional percolation problem. The evaluation of the 
optimal diagram with decreasing values of A is shown in 
Figffl Here one can see that with decreasing A the clus- 
ters become more compact. More precisely one can define 
a finite length-scale, the so called breaking-up length, lb, 
which is rapidly increasing with decreasing A. For small 
A the breaking-up length has been calculated in Ref^ 2 -: 



l b « lo exp 



.4 



(5) 



which is divergent for A — > 0, i.e. in the non-random sys- 
tem limit. In a numerical calculation on a finite sample 



FIG. 1: (Color online) Different type of clusters in the opti- 
mal diagram of the RBPM: spanning clusters [black], clusters 
which touch only the upper (lower) boundary of the strip [red 
(green)] and clusters which have no common points with the 
boundary [yellow]. Lower panel, A = 1/2, standard bond 
percolation; middle panel, A = 5/12; upper panel, A = 4/12. 
Note that the breaking-up length is increasing with decreasing 
A. 



FIG. 2: (Color online) Normalized density profiles, pb{l/L), 
of the RBPM for different widths. The dashed line indicates 
the conformal result in Eq.([6| with the conjectured exponents 
in Eq.([l} and with the boundary parameter a — in Eq.©. 
In the inset the ratio of simulation to theoretical results are 
presented for L — 256 and for two different boundary param- 
eters: a = and a = — 1.0, see Eq.©. 



B. Densities of critical clusters 



of linear size, L, one should have the relation, L 3> lb, 
thus A should be not too small. On the other hand one 
should also be sufficiently far from the percolation limit, 
A = 1/2, in order to get rid of cross-over effects. This 
means that the optimal choice of A is a result of a com- 
promise, which in our case seems to be around A = 5/12, 
when the typical brcaking-up length is about If, ~ 14. 
Most of our studies are made for this value, but in order 
to check universality, i.e. disorder independence of the re- 
sults we have made also a few calculations for A = 21/48, 
too. 



Calculation of the optimal diagram for a given realiza- 
tion of disorder is a non-trivial optimization problem, for 
which very efficient combinatorial optimization algorithm 
have been developed^, which works in strongly polyno- 
mial time. Application of this method made it possible to 
obtain the exact optimal diagram for comperatively large 
finite systems. In order to have an effective strip geom- 
etry we have considered lattices of rectangle shape with 
an aspect ratio of four. The strips have open boundaries 
along the long direction and periodic boundary condi- 
tion was used in the other direction. We mention that 
the same geometry has been used before for percolation, 
too^. The width of the lattices we considered are from 
L = 32 up to 256. Typically we have considered sev- 
eral thousand samples, for the largest system we have 
thousand samples. 



We start to study the density of crossing clusters, 
Pb(l/L), which is given by the probability that a point 
in the position, I, measured perpendicular to the strip, 
belongs to a cluster which touches both boundaries of 
the strip. For percolation in the continuum limit, I 3> 1, 
L ^> 1 and y = l/L the density, pb(y), is calculated 
trough conformal invariancei^: 



x \( Try\ x 
Pb{y) oc (sin Try) |^cos — j 



sin — — — 1 



(6) 

in which Xf, = 5/48 and x s = 1/3 are the scaling dimen- 
sions for percolation^. 

For the RBPM the numerically calculated normalized 
densities, pb(l/L), for different widths are shown in Fig|2] 
All the data fit to the same curve and the finite breaking- 
up length, If,, seems to have only a small effect. 

In the surface region, 1 « i, but I > If, one expects 
from scaling theory: Pb(l) ~ i Is_Ii , which is in accor- 
dance with the limiting behavior of the conformal predic- 
tion in Eq.([6]) with the conjectured scaling exponents for 
Xb and x s in Eq.|T|). In Fig[3]wc have presented pb(l) in a 
log-log plot in the surface region for the largest finite sys- 
tem. Indeed, for I > lb the points are well on a straight 
line the slope of which is compatible with the conjec- 
tured value: x s — Xb = 0.309. We have also estimated 
the asymptotic slope of the curve by drawing a straight 
line through the points in a window [lb + 1/2, lb + 31/2] 
by least square fit. Fixing lb = 15 the estimates with 
varying I seem to have a ~ I 2 correction (see the inset of 
FigEl) and the extrapolated slope is x s — Xb — 0.303(8) 
in agreement with the conjectured values in Eq.|T]). 
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FIG. 3: (Color online) Density profile, p b (l), of the RBPM 
for L = 256 close to the surface in log-log plot. The straight 
(red) line has a slope x s — x\, = 0.309. Inset: estimates of the 
slope using different windows of the fit, see the text. Here the 
full (red) line indicates a parabolic fit. 

We have also checked if the conformal result for perco- 
lation in Eq.© using the conjectured scaling exponents 
for xi, and x s in Eq. © can be used to fit the scaling curve 
for the RBPM for the whole profile. As seen in Figj2]the 
agreement between the numerical results and the formula 
in Eq.© is indeed very good. We have also calculated 
the ratio of the simulation to the theoretical results. In 
this case for the theoretical curve we used the variable: 

y=(l + a)/(L + 2a), (7) 

where a = 0(1) measures the effective position of the 
boundary in the lattice model. By varying a one can 
obtain a better fit in the boundary region, as seen in the 
inset of FigH 

Next we consider the density of those clusters which are 
touching one boundary of the strip, say at y = l/L — > 
0, irrespective to the other. This density, denoted by 
po(l/L), in the continuum approximation is calculated 
for percolation by conformal methods^ as: 

po(y) oc (sin Try) Xh [cos, — J (8) 

This density is analogous to the order parameter profile 
in the system with fixed-free boundary condition a 19 ' 20 . 
The numerically calculated densities are shown in Fig(3] 
for different widths, where we used such a normalization 
that the curves have the same asymptotics at the free 
boundary, i.e. around y = 1. Close to the free boundary 
the densities for different L fall to the same curve, which 
is well described by the conformal formula in Eq.© in 
which we have used the conjectured exponents in Eq.©. 
In larger distance from the free surface the numerically 
calculated densities start to deviate from the conformal 
result for some y < t)l and j/l is a decreasing function 
of L. If we extrapolate the simulated profiles with L the 
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FIG. 4: (Color online) Density profiles, p (l/L), of the RBPM 
for different widths by using such normalization that the 
curves have the same asymptotics around y = 1. The ex- 
trapolated values for l/L > 0.5 are denoted by red squares. 
The dashed line indicates the conformal result in Eq.© with 
the conjectured exponents in Eq.fTJ. Inset: density profiles 
pc.il/ V) for different widths. 



region of agreement with the conformal result is extended 
to the interval 0.2 < y < 1., as seen in FigfJ] 

The finite-size dependence of the densities in this case 
can be attributed to the effect of the finite breaking-up 
length, 4- Close to the touched surface in the continuum 
limit, l b <C I <C L, the density is described by the scaling 
result by Fisher and de Gennes^: po(l) ~ l~ Xb . However 
by approaching the breaking-up length, lb, the increase 
of the profile is stopped and for I < lb po(0 start to 
decrease. This is due to the structure of the connected 
clusters close to the surface. As seen in Fig[T]the number 
of touching sites in a cluster is compcratively smaller for 
the RBPM with A < 1/2 (upper and middle panel of 
Fig© > than for percolation with A = 1/2 (lower panel of 
Fig©. Also for finite widths the small and medium size 
touching clusters are rarely represented for the RBPM. 
By approaching with I the other, free side of the strip the 
crossing clusters start to bring the dominant contribution 
to the density, po(l/L), which is then well described by 
the conformal formula. 

Finally we consider p e (l/L), which is the density of 
points in such clusters which are touching either the 
boundary at / = 1 or at I = L or both. For percolation 
this density is predicted through conformal invariancc 
asii: 

p e (y) cc (smny)~ Xh (9) 

and it is analogous to the order parameter profile with 
parallel fixed spin boundary conditions^. Note that 
we have the relation: p b (y) = p (y) + p\{y) - p e {y), 
with pi(y) — po(l - y). For the RBPM this density 
is strongly perturbed by the finite breaking-up length at 
both boundaries as can be seen in the inset of FigfJ] In 
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FIG. 5: (Color online) Density profile along a vertical line 
with two touching boundary points, p^ ne (l/L), for the RBPM. 
The dashed line indicates the conformal result in Eq.Q with 
the conjectured exponents in Eq.Q. In the inset the same 
quantity is shown for percolation. Here in the analytical for- 
mula in Eq.lO we use xt = 5/48 and x s = 1/3. In both 
figures the boundary parameter in Eq.© is a = 0. 



III. FINAL CLUSTERS IN THE RTFIC 
A. Model 

The random transverse-field Ising chain is defined by 
the Hamiltonian: 



H 



h t a* 



(10) 



4!) 



in terms of the Pauli matrices, o\' z , at site i. The cou- 
plings, Ji, and the transverse fields, hi, are independent 
and identically distributed random numbers. The criti- 
cal point of the system is located at: [ln/i] av = [In J] av , 
where we use the notation [. . .] av to indicate the average 
over quenched disorder. 

We note that the RTFIC is the Hamiltonian vcrsio: 
of the McCoy- Wu model^ 3 ., which is a 2d Ising model 
with layered disorder. In the i-th layer of this model 
the couplings in the vertical and horizontal directions 
are given by: K\(i) and K 2 (i), respectively, which are 
related to the parameters of the RTFIC as: hi = 
— t -1 tanh -1 exp(— 2K\{i)) and Ji = —r~ l Ki(i), where 
in the Hamiltonian limit r — > 0. 



this case we did not try to perform an extrapolation and 
conclude that even larger finite systems would be nec- 
essary to test the conformal predictions in a direct cal- 
culation. In order to try to test the result in Eq.® we 
studied another density which is defined on crossing clus- 
ters, so that one expects to be represented correctly in 
smaller systems, too. Here we define a density, /9g ne (Z/L), 
in crossing clusters and consider points only in such ver- 
tical lines, where at both ends of the given line the clus- 
ter touches the boundaries. Since p^ nc (l/L) is related to 
the operator profile with fixed-fixed boundary conditions 
we expect that it has the same scaling form as the pre- 
viously defined density, p e (l/L). In Fig O we show the 
calculated densities for the RBPM, which is compared 
with the analytical prediction in Eq.®. A similar anal- 
ysis for percolation is shown in the inset of FigO In both 
cases we found that the numerical and analytical results 
for this type of profile are in satisfactory agreement, al- 
though the statistics of the numerical data is somewhat 
low, since just a fraction of ~ L~ 2Xs ~ L" 1 lines can be 
used in this analysis. 

We can thus conclude that all the critical densities we 
considered for the RBPM arc found in agreement with 
the theoretical prediction, which is obtained from the 
corresponding conformal results for percolation by re- 
placing the scaling dimensions with the appropriate (con- 
jectured) values for the RBPM. From an analysis of the 
profile, Pb(l), close to the boundary we have obtained new 
accurate estimate of the critical exponent, x s — x&, giving 
further support of the conjecture in Eq.|T|). 



B. SDRG method 

The RTFIC can be efficiently studied within the frame 
of a renormalization group approac h 35 ' 50 , which is ex- 
pected to lead to asymptotically exact results^. The 
basic feature of this procedure is to successively elim- 
inate those degrees of freedom which have the largest 
local energy scale and thus represents the fastest local 
mode. At a given step of the renormalization the global 
energy scale is defined by f2 = max{ Ji, hi} and the local 
term of value f2 is eliminated from the Hamiltonian. Here 
we have two different elementary renormalization steps: 
cluster formation and cluster decimation. 

i) Cluster formation: if the largest local parameter is a 
coupling, say Ji 3> hi, /13 {hi and /13 being the transverse 
fields acting at the two ends of Ji) then a new spin cluster 
is formed in an effective transverse field: h 23 « hih^/Ji, 
which is calculated in second-order perturbational cal- 
culation. The moment of the new cluster is given by 
M23 = t l 2 + M3j m terms of the moments of the original 
clusters, p,i and /Z3. In the starting Hamiltonian all spins 
have the same moment of unity. 

ii) Cluster decimation: if the largest local parameter 
is a transverse field, say hi ^> Ji,Jz {Ji and J3 being 
the couplings which are connected to the site with hi) 
then the spin cluster is decimated out and an effective 
coupling J23 ~ JiJz/hi is formed between the remaining 
sites. If the decimated spin is at the boundary of an open 
chain no new couplings are formed. 

During renormalization we repeat the elementary dec- 
imation steps, which at the starting period are only ap- 
proximative, but as the energy scale is reduced and the 
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FIG. 6: (Color online) Examples of final clusters in the RTFIC 
for L = 64. Sites in the final cluster are denoted by black 
squares. Upper panel: final cluster in general position; middle 
panel: final cluster containing the boundary spin at £ = 1; 
lower panel: final cluster containing both boundary spins at 
I = 1 and at I = L. 

fixed point, Q* = 0, is approached they become asymp- 
totically exact. In this limit the renormalization group 
equations can be solved analytically. The length-scale of 
the clusters (and bonds) , defined by the linear size of the 
original region which is renormalized to the new variable 
is shown^i to scale as 

£~\n(n/n f , (n) 

where fio is a reference energy scale. On the other hand 
the average cluster moment behaves as^i: 

Ai ~in(n/n )*, <& = ^^. (12) 

Note that the average magnetization at the critical point 
behaves as: m(£) ~ p/£ ~ l~ Xb as lengths are rescaled by 
a factor, £, and Xf, = 1 — <J>/2 is just the scaling dimension 
introduced in Eq.©. 

C. Densities of critical clusters 

Having a finite chain with length, L, we perform the 
decimation until the final cluster, which have a moment 
of fj,(L) -i*/ 2 . Sites of the original chain which belong 
to the final cluster are very strongly correlated and we 
can ask questions about the density of sites in the final 
cluster, i.e. about the probability that a given site is 
contained in a final cluster. The structure of spins in the 
final clusters are illustrated in FigJFo] Note that the final 
cluster in the Id space is disconnected, and the correla- 
tions are realized along the (imaginary) time direction. 
Densities of critical clusters in the RTFIC is studied nu- 
merically. We have considered a large number of (3 x 10 7 ) 
chains of length L = 2 13 = 8192 with open boundary 
conditions. We used the same type of uniform disorder: 
p(u) = 1, for < u < 1 and p(u) = 0, for u > 1, both 
for the couplings and the transverse fields, in this way 
we have satisfied the criticality condition. The strong 
disorder renormalization procedure is performed for each 
chain up to the final spin cluster and then the statistics of 
the sites belonging to the final clusters are investigated. 
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FIG. 7: (Color online) Density profile of the RTFIC consid- 
ering all final clusters of the SDRG procedure. The dashed 
(red) line indicates the formula in Eq.©. Insets: a) Density 
profile close to the surface in log-log plot. The straight (red) 
line has a slope x s — Xb = 0.309. b) Ratio of the numerical 
results and the formula in Eq.®. The parameter in Eq. (JTJ) is 
a=0. (full or red line) and a=-l. (dashed or green line). 



We have studied the density of three different class of 
clusters, which have somewhat analogous definitions to 
the clusters studied for the RBPM. In terms of all final 
clusters (see the upper panel of Fig© we define p{l/L). If 
we consider those final clusters which have the boundary 
point I = 1 (see the middle panel of Fig© we obtain 
po(l/L). Finally, if the clusters contain both boundary 
points, I = 1 and I = L (see the lower panel of FiglfJ]) 
we define pm(l/L). We note that for these densities no 
analytical conjecture is available, since the system is not 
conformally invariant. 

The density of all final clusters is shown in FigJT] 
First we note that close to the boundaries the behav- 
ior of the profile is predicted by scaling theory as p(y) ~ 
y x s -x b> y 4^1, or p(£) ~ P'-**, £ < L. This rela- 
tion is indeed satisfied as shown in inset a) of Fig[7J From 
this inset we can notice that the microscopic length-scale 
of the model, l m , is just a few lattice spacing and for 
I > l m the calculated profile is well described by the 
asymptotic scaling result. The scaling result is valid for 
the formula in Eq. © with the appropriate scaling dimen- 
sions, therefore we tried to compare it with the numerical 
results. As seen in FigJT] the agreement is very good for 
all values of y. To have a more precise check in inset 
a) of FigUJ we have presented the ratio of the numerical 
results and the formula in Eq.©. Here one can notice 
small deviations from unity, which are of the order of 
1%. Consequently the formula in Eq.© is a very good 
fit, however presumably it is not exact. 

We can conclude that the density profiles of the final 
clusters of the SDRG procedure for the RTFIC are well 
described by scaling predictions close to the boundaries. 
The full profiles arc also well approximated by analytical 
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FIG. 8: (Color online) Density profile of the RTFIC consider- 
ing those final clusters of the SDRG procedure which contain 
the site at I = 1. The dashed line indicates the formula in 
Eq. (|13[l . Inset: density profiles at the two surfaces in log-log 
plot. The straight lines have slopes: —Xb = —0.191 (upper or 
green) and x s — x b = 0.309 (lower or red). 



FIG. 9: (Color online) Density profile of the RTFIC consider- 
ing those final clusters of the SDRG procedure which contain 
both boundary sites. The dashed line indicates the formula in 
Eq.@. Inset: ratio of the numerical results and the formula 
in Eq.©. Typical errorbars are also indicated. 



formulae, which are however not fully perfect. 

The density of final clusters which contain the bound- 
ary site at I = 1 is shown in Fig [51 From scaling 
theory one knows the behavior of the profile close to 
the boundaries: po{y) ~ {y)~ Xb , V -C 1 (or p (£) ~ 
(£)-*>, £<^L)andpo(y)^(l-y) x -- Xb , 1 - y < 1 
(or p (L -£)~(L- i) x °- x \ L-£<^L), respectively. 
This behavior is indeed found in the numerically calcu- 
lated profile as seen in the inset of FigJS] The asymptotics 
mentioned above is valid for the formula in Eq.©. We 
tried to fit the numerical results with this formula (with 
the appropriate scaling dimensions), however the weight 
of the the tail at y ~ 1 given by this formula is too large, 
by about a factor of 2. Much better agreement with the 
data can be obtained with the formula: 



Pa(y) oc (sin ?ry) x " (cos ^ 



sin^) +1 



(13) 

which is just the average of the density of clusters which 
touch one boundary an may and may not touch the other 
boundary. As seen in FigJSJ the analytical and numerical 
results arc close to each other for all y, although the 
agreement is certainly not perfect. 

Finally wc consider those final clusters that touch both 
boundaries. The corresponding density, poi(y), is similar 
to the order parameter profile with fixed-fixed boundary 
condition and its functional form for percolation is given 
in Eq.©. The numerically calculated profile is given 
in FigUl Here the comperatively large fluctuations of 
the data points are due to the fact that only a fraction 
~ L~ 2Xa ~ L^ 1 of the samples have a final cluster which 
touches both boundaries. We have compared the calcu- 
lated profile with the analytical formula in Eq.© using 
Xb from Eq.©. The agreement is generally very good, 



but not fully perfect. Small deviations of the order of a 
few percents can be observed (see in the inset of Fig© . 



IV. DISCUSSION 

In this paper we have studied the density of critical 
clusters in two models the critical properties of which are 
dominated by disorder effects. One of the models is the 
two-dimensional random bond Potts model and we con- 
sidered the FK clusters in the largc-g limit. This model is 
expected to be conformally invariant, which means that 
average quantities which are related to FK clusters (such 
as correlation function and magnetization densities) are 
invariant under conformal transformations. We note that 
the RBPM and conventional bond percolation represent 
two different fixed points of the same phase diagram, 
which correspond to < A < 1/2 and A = 1/2 for 
the binary disorder, respectively, see Fig[T] In contrast 
to percolation in the RBPM there is a finite length-scale, 
the breaking-up length, Z^, and results of the continuum 
approximation are expected to hold for lengths which are 
larger than ij. In the strip geometry we have calculated 
the density of points of different type of clusters (cross- 
ing clusters, clusters which touch one boundary of the 
strip, etc.) in analogy with a related study of percola- 
tion ini^. The densities close to free surfaces are well 
described by scaling predictions and from this analysis 
accurate estimate of the critical exponent X g X IS ob- 
tained in agreement with the conjecture in Eq.©. The 
full profiles are compared with analytical formulae which 
are obtained from the corresponding conformal results for 
percolation by using the appropriate values of the bulk 
and surface scaling exponents in Eq.©. We have ob- 
served that the numerically calculated profiles agree well 
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with the conformal results outside the surface region of 
width ~ Ifr. 

The second model we considered is the RTFIC, the 
fixed point of which is expected to control the critical 
behavior of a large class of 2d classical systems with 
anisotropic randomness. Examples are the Ising model 
and the (directed) percolation with layered disorder. In 
these systems scaling at the critical point is strongly 
anisotropic, therefore these systems are not conformally 
invariant. In the RTFIC critical clusters are defined 
through the strong disorder RG procedure. Here spins 
in the final cluster are strongly correlated and play anal- 
ogous role as clusters in percolation or FK clusters in the 
Potts model. The density of final clusters of the RTFIC 
close to the surfaces of the strip are shown to obey scaling 
relations. We also tried to find analytical formulae which 
correctly approximate the numerical profiles. These for- 
mulae, which are borrowed from similar studies of con- 
formal systems have an overall very good description, 



however these are not fully perfect. We have noticed a 
discrepancy of the order of a few percent. 

Our investigations can be extended into different direc- 
tions. For 2d classical systems one can study the density 
of FK clusters in the g-state Potts model, both without 
disorder (for q < 4) and in the presence of disorder (for 
general value of q) . One can also study the density of geo- 
metrical clusters in the 2d random-field Ising mode l 51 ' 52 . 
For the random transverse-field Ising model one possibil- 
ity is to investigate the distribution of final clusters in a 
2d strip. 
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search and Technology under Grant No. ASEP1111, 
by German-Hungarian exchange programs (DAAD-MOB 
and DFG-MTA) and by the Hungarian National Re- 
search Fund under grant No OTKA T048721, K62588, 
M045596. M.K. thanks the Ministere Francais des Af- 
faires Etrangeres for a research grant. 
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